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Decoherence is a major obstacle to any practical implementation of quantum information pro- 
cessing. One of the leading strategies to reduce decoherence is dynamical decoupling — the use 
of an external field to average out the effect of the environment. The decoherence rate under any 
control field can be calculated if the spectrum of the coupling to the environment is known. We 
present a direct measurement of the bath coupling spectrum in an ensemble of optically trapped 
ultracold atoms, by applying a spectrally narrow-band control field. The measured spectrum follows 
a Lorentzian shape at low frequencies, but exhibits non-monotonic features at higher frequencies due 
to the oscillatory motion of the atoms in the trap. These features agree with our analytical mod- 
els and numerical Monte-Carlo simulations of the coUisional bath. From the inferred bath-coupling 
spectrum, we predict the performance of well-known dynamical decoupling sequences: CPMG, UDD 
and CDD. We then apply these sequences in experiment and compare the results to predictions, 
finding good agreement in the weak-coupling limit. Thus, our work establishes experimentally the 
validity of the overlap integral formalism, and is an important step towards the implementation of 
an optimal dynamical decoupling sequence for a given measured bath spectrum. 



I. INTRODUCTION 

In any implementation of quantum information pro- 
cessing by qubits [T] it is crucial to keep the qubits coher- 
ent for long periods of times. The qubits, however, are 
never completely isolated from the environment. This 
coupling to the environment means that after some time 
an entanglement between the system (the qubits) and the 
environment (all other degrees of freedom) is established. 
Usually the environment (bath) consists of many degrees 
of freedom which arc not controlled. This means that one 
has to trace over these degrees of freedom to obtain the 
state of the system. The entanglement to the bath and 
the tracing leave the system in a separable classical state. 
Obviously, the resulting loss of coherence (decoherence) 
is one of the major obstacles towards the successful im- 
plementation of quantum information processing. 

Over the years, several strategies have been suggested 
to cope with the decoherence problem. Clearly, the first 
thing to do is to minimize the coupling to the environ- 
ment. Borrowing ideas from classical error correction 
theory, it was shown that by encoding a logical qubit in 
several physical qubits, it is possible to correct for er- 
rors which are introduced in the calculation process [1]. 
Quantum error correction protocols can correct errors up 
to some maximal rate. The upper bound on the error 
probability of a quantum gate depends on the details of 
the error correction protocol, but the typical values range 
from lO"'^ to IQ-"^. The gate can also be an identity gate 
which describes the storage of information in a quantum 
memory. Such a memory is essential in the architecture 
of a quantum network in order to enable scalability . 

Dynamical decoupling (DD) is a technique that was de- 
veloped to further reduce the error rate below the fault 
tolerant threshold [^-j^- The main idea of DD is to use 
external control fields that induce rotations of the qubit 
in the Bloch sphere such that the overall decoherence is 



reduced. This method was first considered in the con- 
text of nuclear magnetic resonance [51-[TU] . In the field of 
quantum information processing, the pursuit for fault tol- 
erance has pushed forward the development of a similar 
formalism for controlling the decoherence of noisy qubits 
[31 m |S] . The simplest example of DD is the well-known 
Hahn echo technique: a single population inverting pulse 
(•TT-pulse) is introduced exactly at half the final observa- 
tion time [TT]. The echo technique is widely used and 
very efficient in counteracting the effect of a quasi-static 
coupling to the environment. Once the coupling spec- 
trum extends to higher frequencies, the echo technique 
fails and more elaborate control sequences are needed. 

A strategy aimed at maximizing the decoupling has 
been developed by Kofman and Kurizki P^^HTT] , based on 
a formula that relates the decoherence rate to the over- 
lap of the bath coupling spectrum on the control power 
spectrum. This overlap can be minimized under the con- 
straint of the available control field energy [121 H] • In 
order to successfully implement this decoupling control 
strategy, detailed knowledge of the bath-coupling spec- 
trum is required. In this work, we explain how to measure 
this spectrum, and then use it to calculate the coherence 
of the system after a DD sequence has been applied. The 
experiments we present are performed in an ensemble of 
ultra-cold atoms held together by a potential induced by 
a far-off resonance laser field. This system is used not 
only for demonstration purposes, but also because it can 
serve as a quantum memory |18H23j . 

Two of the most widely used physical qubits are pho- 
tons and neutral atoms. Photons are easy to produce, 
manipulate and transport. They interact weakly with 
their environment and therefore can remain coherent for 
long travel distances. This last advantage is also their 
disadvantage: interactions between photons are usually 
very small, making the implementation of an all-optical 
two-qubit gate very difficult. Atoms, on the other hand. 
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are easy to keep in one place, and can interact strongly 
with other atoms and electromagnetic fields. It is there- 
fore sensible to use atoms as "stationary qubits" for stor- 
age and manipulation and use photons as "flying qubits" 
that carry the information between distant sites. 

One of the controlled schemes of interaction between 
atoms and photons is electromagnetically induced trans- 
parency (EIT) [23]. In this scheme, atoms with a lambda- 
shape energy structure interact with two light fields 
called "pump" and "probe" . The pump is usually much 
stronger than the probe, and is used to control the inter- 
action strength between the probe and the atomic ensem- 
ble. Turning off the pump while the probe is propagating 
in the atomic ensemble leads to a conversion of the pho- 
tonic excitation into the coherence between the two low- 
lying states of the atoms. This is sometimes called "stor- 
age of light" , although only the coherence which was car- 
ried by the light is actually stored in the ensemble. The 
beauty of this conversion process is that it is reversible, 
which makes the atomic ensemble a true memory. 

In order to increase the efficiency of the storage and 
retrieval processes, it is desirable to work with atomic 
ensembles with high optical depth [23 [55] . This is be- 
cause the coupling of the atoms to the external electro- 
magnetic field scales as the square root of the number of 
atoms in a volume where the light intensity is approxi- 
mately uniform. Working at high optical depth, however, 
usually implies that the atomic density is high and the 
rate of inter-particle collisions is large compared to the 
storage time. The coherence properties of the atomic 
ensemble are markedly changed due to elastic collisions 
[27j . From the point of view of a particular atom, other 
atoms can be regarded as the bath. The collisions with 
other atoms define the nature of the coupling to the bath. 
The three important physical quantities in the descrip- 
tion of the coUisional bath are: the collision rate, the 
inhomogeneous dephasing rate (i.e. the dephasing with- 
out collisions) and the harmonic confinement oscillation 
frequency. The spectral behavior of the coUisional reser- 
voir stems from the interplay between these quantities. 

The goal of this work is to present a general method of 
measuring the bath-coupling spectrum, and demonstrate 
its usefulness in calculating the result of any control field. 
We start by reviewing the mathematical formalism of 
the overlap integral spectrum and the power spectrum 
of the control sequence that acts as filter function (Sec- 
tion |ll| . Using this formalism it is possible to calculate 
from the bath-coupling spectrum the coherence at some 
observation time with any control sequence, as long as 
the so-called 'weak coupling' limit applies. In particular, 
it can be used with a null control sequence, which simply 
describes a dephasing process. For a known filter func- 
tion of the control field, it is possible to de-convolve the 
overlap integral to obtain the bath-coupling spectrum. 
We present this technique for a constant-power contin- 
uous control field in Section [HT] and illustrate it by full 
3D Monte-Carlo simulations. We then employ the tech- 
nique to measure directly the coUisional bath spectrum 



in our cold atomic ensemble (Section IV ) . The measured 



bath shows interesting non-monotonic behavior due to 
the oscillatory motion of the atoms in the trap. Once the 
bath-coupling spectrum is known, it is interesting to test 
its implications for the outcome of other DD sequences. 
This is done in Section ]Vj where we use the measured 
bath and the overlap-integral formalism to calculate the 
coherence under the well-known CPMG, UDD and CDD 
sequences. We then experimentally apply these DD se- 
quences and compare the results with the predictions of 
the theoretical calculations based on the measured bath. 
We conclude and give our outlook in Section [yl] 



II. THE SPECTRAL OVERLAP INTEGRAL 
APPROACH TO THE FIDELITY CALCULATION 



In this section we review the formalism [T2H17] which 
is useful to calculate the ensemble coherence at a given 
time. Here we give a simplified version of this approach 
assuming an effective two-level stochastic Hamiltonian 
weakly coupled to the environment. 

We consider the same model as was described in Ref. 
|28) . Our ensemble consists of two-level systems (TLS), 
and due to inhomogeneities in the external environment, 
the transition energy of each TLS is different than the 
frequency it may have in free space (we shall call this 
difference the detuning). We reduce the full many-body 
Hamiltonian to an effective single particle Hamiltonian 
given by 

H ^h[ujo + 5{t)]\2) {2\ + hVL{t)\2) {l\ + h.c. , (1) 

where wq is the free space transition frequency between 
the two states, 5{t) is the detuning, and Vl{t) is a classical 
external control field which is going to be used for the 
DD. The detuning is a random function of time, and we 
assume it is averaging to over different realizations of 
the Hamiltonian. The reduced density matrix can be 
calculated by 



p{t) = m)) m)\ 



(2) 



where 4){t) is the state of the system at time t with a spe- 
cific realization of the stochastic Hamiltonian, and the 
average is taken over many realizations. \4>{^)) is the 
initial state of the system, and ultimately our goal is to 
conserve this state in the ensemble for as long as possi- 
ble. To characterize how well this goal is achieved it is 
instructive to introduce the fidelity function: 



F = {m\p{t)\m) , 



(3) 



which starts at 1 and decays to 0.5 at long times. 

To make the model clearer, we now relate these defini- 
tions to atomic ensembles [27]. In our system, ultracold 
atoms are confined by an external optical potential. Al- 
though the atoms have many energy levels, we confine our 
attention to two low-lying states with negligible sponta- 
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neous decay. The potential is not exactly the same for 
the two states. The main source inhomogeneous broad- 
ening is that the energy difference between the two sates, 
which is proportional to the total energy of the atom, is 
changing in space 29 . The phase space density of our 
ensemble is low enough to be considered classical (far 
from quantum degeneracy). The motion of the atom in 
the external potential causes the transition frequency to 
change in time. Since the ensemble consists of many 
atoms, each with a different trajectory, measuring the 
coherence of the ensemble gives the averages discussed 
above in a single shot. In the experiment, the coherence 
is measured in a time-domain Ramsey-like experiment, 
as described in Ref. [27j . 

We shall define two functions which are going to be 
useful later for the calculation of the fidelity. First, the 
bath coupling spectrum characterizes the spectral con- 
tent of the coupling of the system to the bath and is 
defined as 



G{f) 



-2771 fr 



{S{t) ■ S{t + t)) dT 



(4) 



where (...) stands for the averaging over many realiza- 
tions of d{t). G{f) is the Fourier transform of the time 
correlation function of S(t). The second function we de- 
fine is a 'filter function' which characterizes the power 
spectrum of the control sequence: 



Ftif) 



'2iTift ^ 



cos I / rt{T)dT ) dt 



(5) 



Note that although Ft{f) is defined in the frequency do- 
main, it does depend on the observation time t. 

A simple expression for the fidelity at time t can be 
obtained under the weak coupling assumption, namely 
that the fidelity decay during the bath correlation time is 
negligible. In this case the fidelity is given by the overlap 
integral [HUmi]: 



1 
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'R{t)t 



(6) 



where R{t) is the decay rate given by the overlap integral 



R{t) 



2a 



G{f)Ftif)df 



(7) 



where a is a constant between and 1, depending on 
the assumptions for the statistics of the initial state, in 
particular, a = ^ for an initial state which is an equal 
superposition of the two internal states. Generally speak- 
ing, R{t) is not constant and therefore the decay in not 
exponential: a signature of the non-Markov time-domain. 

A related quantity is the ensemble coherence, which 
is defined to be the normalized off-diagonal element of 
the reduced density matrix |30| . In the overlap integral 
framework, the coherence is given by C{t) = e-«(*)*. In 
the Bloch sphere representation, the coherence is given by 



the length of the Bloch vector, which is usually measured 
by quantum state tomography [1] [28] . 



III. 



HOW TO MEASURE THE BATH 
COUPLING SPECTRUM ? 



In order for Eqs. (6|7) to be useful, one has to know 
the bath coupling spectrum G{f), preferably from ex- 
perimental data. This spectrum can be inferred from 
Eq. (0 by measuring the coherence for a given pulse se- 
quence with a known Ft{J). Obviously, if Ft{f) where a 
Dirac S function, the decay rate would be linearly propor- 
tional to the bath coupling spectrum. A good approxi- 
mation to a (5 function can be obtained if the control field 
is nearly-continuous and on-resonant field. This control 
field causes Rabi oscillations with a frequency /o which 
depends on its strength. The resulting filter function is 
a sine- function centered around /q. To be more precise, 
the filter function is given by 

Ft{f)^^t[sinc\t{f~fo))+smc\t{f + fo))] , (8) 

where t is the pulse duration. For this pulse, assuming 
/qT ^ 1, Eq. ([T]) yields the following decoherence rate: 



(9) 



We see, then, that by scanning the strength of the con- 
trol fields (thus scanning /o) and measuring the decay 
of coherence, it is possible to directly measure the bath 
coupling spectrum. 

To illustrate the method and check its sensitivity to 
the weak-coupling assumption, we have run a 3D Monte- 
Carlo simulations of a cold atomic ensemble trapped in 
an optical potential. The simulation solves for the clas- 
sical Newtonian motion of 3500 atoms. The initial con- 
ditions of the atoms are drawn from a Boltzmann distri- 
bution, assuming a temperature of 7^K. As explained 
previously, the fluctuations in our system originate from 
elastic collisions between the atoms, and it is therefore 
necessary to include them in the simulation. In order not 
to run into computational complexity problems, we have 
developed a mean-field technique [3T]. Once the trajec- 
tories of all the atoms are calculated, we calculate the 
energy shift of the internal states induced by the exter- 
nal potential along the trajectory of each atom. In our 
7 — 1.06/im-wavelength dipole trap, the differential shift 
of the internal states is 6.6 • 10~^ times the overall po- 
tential [29 . These energy shifts are then used to solve 
the Bloch equations [32] in the presence of the control 
field. We calculate the decoherence rate by computing 
the length of the Bloch vector at different times and fit 
it to a decaying exponent. 

We simulate a measurement of the bath spectral func- 
tion as described above by a continuous control field, in 
two parameter regimes which correspond to weak and 
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strong coupling. The results of the simulations are de- 
picted in Figure [l] At low frequencies, the spectrum fol- 
lows a Lorentzian curve which is expected since the fluc- 
tuations originate from a scattering process with Poisson 
statistics 28J. Also apparent in the spectrum is a non- 
monotonic feature, due to the oscillatory motion of the 
atoms in the trap (for more details see the next section). 
The only difference between the two simulations is the 
number of atoms, resulting in different correlation times 
of the bath. It is interesting to note that by increasing 
the number of atoms one also decreases the decoherence 
rate at low frequencies, a result of coUisional narrowing 
(see Ref. [27] for more details). Thus, increasing the 
number of atoms increases both the collision rate and 
the coherence time, thereby pushing the ensemble deeper 
into the weak-coupling regime. Note that the relevant 
decoherence rate for considering whether the system in 
weakly coupled is the decoherence rate obtained with the 
decoupling pulse. This means that a system can be in 
the strong-coupling regime at low frequencies and in the 
weak-coupling regime at high frequencies, whence devi- 
ations from the calculated spectrum are larger at low 
frequencies. 
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IV. MEASUREMENT OF THE BATH 
COUPLING SPECTRUM 

We have used the technique described in the previ- 
ous section to measure the bath coupling spectrum in 
an ensemble of colliding ultra-cold atoms. In the exper- 
iment, "^^Rb atoms are trapped in an optical potential 
created by a far-off-resonance laser with a wavelength of 
1.064/Lim. Initially ~ 10^ atoms are trapped and cooled 
in a magneto-optical trap, and further cooled by the 
Sisyphus [33 and Raman-sideband techniques [34]. We 
then use rapid adiabatic passage with a constant RF ra- 
diation and a ramped-up magnetic field to transfer the 
atoms from the state |5^S'i/2, F = 1; = l) to the state 
|5^S'i/2^ — Imp ~ — l)for ^ 69% of the atoms, and the 
rest to the state |5^S'i/2,-F' — ^'^mf = O). The thermo- 
dynamic parameters of the ensemble are measured by 
absorption-imaging of the cloud. 

The optical trap consists of two beams crossing at an 
angle of 28° after passing through a zoom system ca- 
pable of controlling their waist. We start by collecting 
the atoms in a 180/im-waist beam, and then dynamically 
compress the waist down to 50/im. The polarization of 
the two crossing beams is parallel to the magnetic field, 
and their frequency differs by 120MHz, in order to elim- 
inate standing waves. In this experiment we choose to 
work with the two states |1) ~ |5^5i/2^ — Imp = — l) 
and |2) = |5^5'i/2F = 2mp — At the applied mag- 

netic field of 3.23 Gauss these states are Zeeman in- 
sensitive to magnetic field fluctuations to the first or- 
der [35]. The two states are separated by an energy of 
2tt11 ■ 6.833GHz, but are slightly affected by the differen- 
tial AC Stark shift of the dipole trap [33]. As explained 



FIG. 1. Numerical simulations of the bath-coupling spectrum 
inferred from decoherence measurements as compared to the 
direct calculation. We simulate the classical motional of Rb 
atoms in the trap. The blue solid line is the bath spectrum as 
calculated directly from the detunings along the atomic tra- 
jectories using Eq. Q. The red circles are the decoherence 
rates of atoms subject to a continuous field (long pulse) with 
the Rabi frequency plotted on the x-axis. The conditions of 
the simulations are a temperature of T = TfJ^K, radial and 
axial oscillation frequencies of 2tt ■ 600Hz and 2t: ■ 160Hz, re- 
spectively, and number of atoms, N = 2 ■ 10* for graph (a) 
and A'^ — 1.5 • 10^ for graph (b). The deviation of low fre- 
quencies of the simulated decoherence measurements from the 
directly simulated spectrum is due to the breakdown of the 
weak-coupling assumption. The number of atoms is differ- 
ent in the two graphs and so is the frequency at which the 
breakdown occurs. In graph (a) the non-monotonic spectral 
feature (peak) is smeared due to the much larger Lorentzian 
width. 



above, for a moving atom this shift is time-dependent 
and follows the trajectory. 

Since Am = ±2 between the two internal states, the 
external control il{t) has to be effected by a two-photon 
transition. We employ RF radiation at 2.15MIIz and mi- 
crowave radiation at 6.832527928GIIz, which is chosen 
such that both frequencies arc detuned by 90kHz from 
the intermediate level |5^5i/2-F = ^mp = O). The max- 
imum Rabi frequency we achieve is $7 = 27r • lOOOHz. 
To detect the state of the atoms we use a state-selective 
fluorescence-detection scheme, similar to the one de- 
scribed in Ref. [55] . 

Two main difficulties arise from driving the system by 
a strong control field. First, at high Rabi frequencies any 
inaccuracies from shot to shot in the field strength trans- 
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FIG. 2. Calibration of the microwave dressing effect. The 
y-axis shows the population at F—2 as a function of the 
control pulse frequency shifted by fc, the center frequency 
without the dressing. The data was taken with a Rabi fre- 
quency of SI = 340Hz. We fit the data to the theoretical curve 
r[l + cos(<^ + Ty(i^ 



A 



fo)^ + P)], where T is the 



(f-fo)^+P 

duration of the Rabi pulse and / is the field frequency. The fit 
gives fo, the dressing for the applied Rabi frequency. In this 
example, the shift is — 37Hz. The inset shows the microwave 
dressing shift as a function of the Rabi frequency squared. 
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FIG. 3. Raw data of coupling spectra. The points are the 
population at F=2 as a function of the control pulse Rabi fre- 
quency. The experiment was done for 160k atoms in a trap 
with a radial and axial oscillation frequencies of 27r-910Hz and 
27r-240IIz, respectively. The field dressing results in the reduc- 
tion of the contrast which is dominant in high frequency. On 
the right we show that after correction for this field-dressing 
the contrast does not decay at higher frequencies even after 
Isec. The contrast between the two graphs is different since 
the observation time is 0.4sec for the left graph and Isec for 
the right graph. The difference in the mean population stems 
from m-changing transitions in the F — 2 hyperfine level. The 
coherence is calculated by a maximal likelihood method (see 
text), and is determined by the envelope of the scattered data 
points. 



late into noise in the atomic population. For example, 
an average noise of 1% in the Rabi frequency transforms 
into ^ 100% noise in the atomic population, if we drive 
our system at = 27r • lOOHz for Isec. Strictly speaking, 
this noise in the control field results in a reduction of the 
fidelity. In order to avoid this difficulty, we randomize 
completely the initial population difference by an addi- 
tional pulse that rotates the atomic state by a random 
uniformly distributed angle and analyze ~ 30 data points 
taken for the same pulse power using envelope spec- 
troscopy. The corresponding Bloch-vector length is es- 
timated using the maximum likelihood estimator. Since 
the angle of the Bloch vector is now completely random 
we can assume that its z-component is Csin(<i>), where 
<i> is a uniformly distributed random phase. The results 
of the maximum-likelihood algorithm that extracts the 
Bloch vector length C coincide with the extreme values 
of the measured z-component (the population inversion) . 

The second difficulty is due to the energy-level shift 
induced by the control field [37 . This "field dressing" of 
the levels depends on the microwave field strength (we do 
not observe the shift due to the rf field) . When the bath 
spectrum is measured using a continuous control field, 
this effect causes the control to have a detuning which 
depends on its strength (and therefore on its Rabi fre- 
quency). We measure the shift by applying a long pulse 
which induces Rabi oscillations and detect the popula- 
tion at state |2) while scanning the frequency of the mi- 
crowave field. Such a measurement is shown in Figure [2] 
This data yields the relative-level shift for a given Rabi 
frequency. We then repeat this measurement for several 
Rabi frequencies and obtain a calibration graph, shown 
in the inset of Figure [2] As expected from theory, the 



dressing at small field strengths is linearly proportional 
to the control- field Rabi-frequency squared [37]. 

In Figure [3] two examples of raw data from a bath- 
spectrum measurement are shown. The scattering of the 
data points is a consequence of our randomization tech- 
nique explained above. The bath-coupling spectrum at 
a given frequency is related to the envelope of the scat- 
tered points. The data presented in Figure [3^ was taken 
without changing the control field frequency (not to be 
confused with its strength which controls the Rabi fre- 
quency) . 

The second difficulty described above leads to a sub- 
stantial reduction of the contrast at higher Rabi frequen- 
cies due to the detuning of the control field from reso- 
nance, caused by the strong field-dressing of the levels. 
In order to eliminate this effect, we change the frequency 
of the microwave field as its power changes to compensate 
for the field dressing. The result is shown in Figure [3|d to 
restore contrast at high frequencies. In general, this cor- 
rection should also be applied in other DD sequences. In 
the cases we present below, however, the dressing effect 
is very small and unimportant. 

To obtain the bath-coupling spectrum, we repeat the 
measurement described above with at least three differ- 
ent Rabi pulse durations. From each such experiment we 
get the coherence at that time for all Rabi frequencies. 
For each Rabi frequency we fit the decay of the coherence 
at different times to a decaying exponent, from which we 
obtain the decoherence rate. As explained above, the 
fitting to an exponent is justified in the weak-coupling 
regime. In Figure|4]we plot the experimentally measured 
decoherence rate as a function of frequency. Similarly to 
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FIG. 4. Comparison of the measured bath spectrum (plus 
signs) to the simulated spectrum (soUd Une) and decoherence 
rate (circles). The conditions for the experiment and simu- 
lation are 340fc atoms at a temperature of 3.5^K, and with 
the power of the trapping laser being 1.7W, creating a trap 
with a radial and axial oscillation frequencies of 2n ■ 450Hz 
and 2n ■ 120Hz, respectively. The simulated spectrum was 
shifted by a constant of 0.4s~^ to better fit the data. This 
shift accounts for Ti processes and a bias which arises due to 
noise in the measurements. 



FIG. 5. The coUisional bath-coupling spectrum, the fil- 
ter functions and their overlap, calculated from the formula 
(Eq. Q) for two DD sequences: CPMG and UDD. The dash- 
dotted line is a Lorentzian spectrum which describes the Pois- 
son statistics of elastic collisions. The filter function for the 
CPMG decoupling sequence (solid line) and UDD (dotted 
line), both consisting of 16 vr-pulses for an observation time 
of Isec. The UDD spectral overlap is larger at lower frequen- 
cies compared to the that of CPMG. Since the bath-coupling 
spectrum at these frequencies is higher (which means a faster 
decay), the resulting decoherence rate for UDD is larger than 
for CPMG. 



the simulations presented in Figure[T] the measured spec- 
trum follows a Lorentzian at low frequencies and shows a 
non-monotonic feature at twice the axial oscillation fre- 
quency. This feature arises due to the rapid oscillation 
of the atoms in the potential. The reason is that the 
detuning function S{t) always has a component at twice 
the oscillation frequency of the trap (since the detuning 
is proportional to the potential which scales as a;^). Mea- 
suring at twice this frequency couples to this component 
and produces higher decoherence |38) . 

For comparison we plot in Figure |4] the results of our 
Monte-Carlo simulation with the same parameters as 
those measured in the experiment. The agreement be- 
tween the measured and simulated spectrum is quite sat- 
isfactory. We observe some shift between the two spec- 
tra in the position of the non-monotonic feature. This is 
probably due to small inaccuracies in the measured laser 
power and waist used in the simulation. The measured 
width of the second peak is larger in the experiment, 
probably due to the anharmonicity of the Gaussian con- 
fining potential which is not simulated. 

Figures [T]|4] prove the ability of this technique to ob- 
serve small dynamical effects through the spectrum. We 
note that although other methods to measure the oscil- 
lation frequency of the trapping potential exist, the use 
of the spectrum to infer this quantity is unique in that it 
does require excitations of the atoms. 



V. DYNAMICAL DECOUPLING AND THE 
BATH SPECTRUM 

As explained in Section [nj the bath coupling spectrum 
is most useful for calculating the outcome of any DD se- 
quence. In general, for a given bath spectrum and a set 
of constraints (such as the maximum duration, energy or 
action of the control pulse), it is possible to devise an 
optimal decoupling sequence [17l |39j. There are, how- 
ever, well-known sequences which are suitable for certain 
classes of noise spectra. Here we concentrate on three 
such sequences: the Carr-Purcell-Meiboom-Gill sequence 
(CPMG) [lOHSg, the concatenated DD sequence (CDD) 
[B], and the sequence developed by Uhrig (UDD) Al- 
though in this work we restrict ourselves to sequences 
of brief 7r-pulses, this is not required by the theory. To 
compare the performance of the three sequences we con- 
sider the decoherence they yield for a given number of tt 
pulses. 

In Figure [s] we depict schematically the Lorentzian 
spectrum of our collisional bath and the filter function 
(Eqs. ^-^) of the CPMG and UDD sequences with 
16 TT pulses. The CPMG spectrum has a sharp peak at 
/o = SHz, half the pulse rate, and smaller peaks at lower 
and higher frequencies whose envelope is proportional to 
!/(/ ~ /o)^- There are also peaks at harmonies centered 
around 3/o but their area is 9-times smaller than that of 
the main peak. UDD has exponentially reduced coupling 
at low frequencies. This can be very useful for bath spec- 
tra which have a cutoff at high frequencies [43j . which 
is not the case of the Lorentzian spectrum. Namely, 
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the UDD exponential suppression of the lower sidebands, 
comes with a price: for the same number of pulses the 
maximum coupling is at lower frequency, hence for many 
physical coupling spectra, including Lorentzian spectra, 
it will sample portions of the spectrum with larger deco- 
herence. CDD (not shown in the graph) has similar prop- 
erties to UDD. It has an exponentially reduced coupling 
in part of the spectrum, but its overlap with a Lorentzian 
is overall the same. 

Experimentally, we apply the DD sequences for a du- 
ration of 0.4sec, on an ensemble prepared with the same 
conditions as in the previous section. The atoms are ini- 
tialized to an equal superposition \'^) — |1) -t- |2), which 
we want to preserve. The duration of a 7r-pulse was cho- 
sen in these experiments to be 2.3msec, corresponding to 
a maximal Rabi frequency of 217Hz. To reduce the noise 
introduced by the inaccuracies in the control field, we 
switch the phase of the control by a 7r-phase every two 
pulses TT, TT, — TT, — TT, ... |28] . The final atomic state at the 
end of the sequence is determined using state tomogra- 
phy. Our measurement of the population gives only the 
z-component of the Bloch vector. In order to measure the 
other two components we apply before the measurement 
an additional 7r/2 rotation around the x oi y axes. We 
tune the duration of our tt pulses in the DD sequences 
such that the there is only a small z-component after the 
sequence. The x and y components are measured by a 
Ramsey technique, in which we vary the phase of a fi- 
nal 7r/2 pulse and measure the contrast of the fringe. In 
other words, the state tomography gives us the ability 
to tune correctly the vr pulse duration such that the fi- 
nal state always resides in the equatorial plane of the 
Bloch sphere. It is then more accurate to measure the 
coherence by applying a final 7r/2 pulse with a controlled 
phase and detect the population in the z-axis. Finally, we 
calculate the coherence time by fitting the decay of the 
fringe contrast at different times to a decaying exponent. 

The results of these measurements are presented in 
Figure [6| For all DD sequences the coherence time in- 
creases as the number of pulses grows. As expected from 
the overlap between the filter functions of the DD se- 
quences and the bath spectrum (see Figure[5]) , the CPMG 
surpasses both UDD and CDD, which in turn behave 
quite similarly. We stress that our ranking between the 
DD sequences is only valid for our experimental colli- 
sional bath. Even more interesting than the relative per- 
formance of the DD sequences is to compare them with 
the theory presented in Section |ll] To this end, we take 
the measured bath spectrum from the previous section 
(Figure [4| and use it in the overlap integral formalism to 
calculate the predicted outcome of each sequence. There 
is, however, a depopulation Ti process in our system 
which results from m-changing transitions in the F=2 
hyperfine level. We measure this timescale directly to be 
Ti = 2.2sec in our experimental conditions. Since the Ti 
and T2 processes are not correlated, the total decay is 
the sum of the decay rates of both processes. In practice, 
we subtract from the measured bath spectrum the bias 
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FIG. 6. Coherence times under different DD pulse sequences. 
We compare the sequences CPMG (blue squares), UDD 
(green triangles) and CDD (black diamonds) for a given num- 
ber of pulses at an observation time of 0.4sec. The perfor- 
mance of both UDD and CDD is worse than the CPMG, as 
expected from their overlap with the bath spectrum (see Fig- 
ure |5|. The solid lines is the coherence time as calculated 
based on the overlap integral of the sequences filter functions 
with the measured bath spectrum. The conditions of the en- 
semble are the same as for the measured bath spectrum in the 
previous section. Note that the CDD sequence cannot have 
any number of pulses, and this is the reason why there are 
only 5 data points for this sequence. 



which results from the Ti process and the noise in the 
measurement. By doing so we can calculate only the de- 
coherence resulting due to the fluctuations of 5{t) (the T2 
process decoherence rate). We then add to it the decay 
due to the Ti process, which is given by 2/Ti. The total 
decoherence rate is depicted as a solid line in Figure |6j 
and agrees very well with the coherence times that were 
measured directly. We attribute the small deviations at 
low frequencies to a departure from the weak-coupling 
regime. 



VI. CONCLUSIONS AND OUTLOOK 

To summarize, we have presented a method of measur- 
ing directly the bath-coupling spectrum, and used it for 
optically trapped ultracold atoms. The measured spec- 
trum follows a Lorentzian shape at low frequencies, which 
is expected, since the source of fluctuations is s-wave scat- 
tering between the atoms [21] ■ At higher frequencies, 
though, the spectrum exhibits non-monotonic features 
which arise due to the rapid oscillatory motion of the 
atoms in the trap. The usefulness of the concept of the 
bath-coupling spectrum is tested by its ability to pre- 
dict the outcome of any DD control sequence. We have 
measured directly the performance of three well-known 
sequences: CPMG, UDD and CDD. The comparison of 
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these measurements to a calculation based on the overlap 
integral with the directly measured bath spectrum proves 
the validity of this framework at the weak-coupling limit. 
As was already shown numerically in [2HI , the CPMG se- 
quence is found to be superior to the UDD and CDD 
sequences, for our coUisional bath spectrum. 

As noted above, the weak-coupling assumption is es- 
sential for the overlap-integral formalism to be correct. 
Under certain experimental conditions and at low fre- 
quencies, our bath does not fulfil this assumption, and 
indeed we observe deviations from the corresponding the- 
oretical calculations. Measuring the coupling spectrum 
under these conditions is a major challenge. One way 
of addressing it is by using a sequence that mostly re- 
sides in the weak-coupling part of the spectrum, while 
a small sideband probes the low frequencies of the spec- 
trum. The idea is to keep most of the overlap between 
the spectrum and the filter function in the weak-coupling 
domain, while still probing the strong coupling domain 
of the spectrum. An example of such a pulse is: 

nit) = 27r/o(l + l3^cosi27rUt)), (10) 
Jo 

where /o is the carrier Rabi frequency, /3 is the modu- 



lation index and fm is the AM modulation frequency. 
By measuring the decoherence twice: once with the side- 
band and once without, it is possible to extract the cou- 
pling spectrum at the sideband frequency, even if this 
frequency is in the strong-coupling part of the spectrum. 

For a given bath spectrum and a set of constraints on 
the control field, an optimal decoupling sequence can be 
constructed [T71 • This is done by solving the Euler- 
Lagrange equation which minimizes the overlap of the 
control-field filter function and the bath-coupling spec- 
trum. This procedure is yet to be demonstrated exper- 
imentally. In our system, as well as in other systems, 
the spectrum is non-monotonic and as such we expect a 
non-trivial optimal decoupling sequence. What can com- 
plicate this picture is noise in the control field. In this 
work we have circumvented this issue by using an enve- 
lope spectroscopy method. It is possible to extend the 
overlap-integral formalism to include the classical noisy 
control as a second (classical) bath-like spectral function. 
The decay rate is then given by two overlap integrals 
which can be solved to find the optimal DD sequence 
for a noisy control. Since noise in the control is always 
present, this extension of the theory is both necessary 
and practical. 

We acknowledge the financial support of MIDAS, MIN- 
ERVA, GIF, ISF, and DIR 
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